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ABSTRACT 

The detection of diffuse radio emission associated with clusters of galaxies indicates populations of rela- 
tivistic leptons infusing the intracluster medium. Those electrons and positrons are either injected into and 
accelerated directly in the intracluster medium, or produced as secondary pairs by cosmic-ray ions scattering 
on ambient protons. Radiation mechanisms involving the energetic leptons together with decay of neutral pions 
produced by hadronic interactions have the potential to produce abundant GeV photons. Here, we report on 
the search for GeV emission from clusters of galaxies using data collected by the Large Area Telescope (LAT) 
on the Fermi Gamma-ray Space Telescope {Fermi) from August 2008 to February 2010. Thirty-three galaxy 
clusters have been selected according to their proximity and high mass, X-ray flux and temperature, and indica- 
tions of non-thermal activity for this study. We report upper limits on the photon flux in the range 0.2-100 GeV 
towards a sample of observed clusters (typical values 1-5 x 10"^ ph cm"^ s~') considering both point-like and 
spatially resolved models for the high-energy emission, and discuss how these results constrain the character- 
istics of energetic leptons and hadrons, and magnetic fields in the intracluster medium. The volume-averaged 
relativistic-hadron-to-thermal energy density ratio is found to be < 5-10 % in several clusters. 

Subject headings: cosmic rays — galaxies: clusters: general — gamma rays: galaxies: clusters — radiation 
mechanisms: non-thermal 
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1. INTRODUCTION 

Clusters of galaxies are the largest viriaUzed structures in 
the Universe. In addition to a high concentration of galaxies, 
they contain ionized gas distributed in the intracluster medium 
(ICM) which emits thermal bremsstrahlung radiation in the 
soft X-ray (SXR; 2-10 keV) energy range. The dynamics of 
galaxies and thermal gas, together with gravitational lensing 
studies, indicate a dynamically dominant dark matter com- 
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ponent. Finally, many clusters show a 'non-thermal' (NT) 
emission (e.g. radio synchotron) component signifying pop- 
ulations of relativistic leptons in combination with an appre- 
ciable magnetic field, and necessarily, turbulence and shocks 
permeating the ICM. 

Galaxy clusters are also unique reservoirs of cosmic ray 
(CR) hadrons. Sparse thermal gas (n ^ lO"-'- 10"^ cm"^) al- 
lows energetic ions to survive over cosmological timescales, 
leading in principle to a situation in which an appreciable frac- 
tion of the total cluster pressure is provided by accumulated 
CRs confined to the ICM by turbulence or chaotic magnetic 
fields (Volk et al. 1996; Berenzinsky et al. 1997). This sce- 
nario can be directly tested by gamma-ray telescopes (e.g. 
EnBhn et al. 1997; Colafrancesco & Blasi 1998; Berring- 
ton & Dermer 2004; Pfrommer & EnBlin 2004) because high- 
energy photons are expected to be produced either through in- 
elastic proton-proton collisions and subsequent neutral-pion 
decay, or by processes involving energetic electron-positron 
pairs produced as secondaries of the hadronic interactions. 
CR leptons can also be injected by galaxies and active galac- 
tic nuclei and re-accelerated by turbulence and/or merger or 
accretion shocks in the ICM. In contrast to the CR ions, pop- 
ulations of relativisitic leptons are expected to be shorter-lived 
(< Gyr; Petrosian 2001) and depend more upon recent merger 
events than on the full formation histories of their host clusters 
(Blasi et al. 2007 and references therein). 

The presence of NT activity in the ICM of some clusters 
was first discovered as diffuse (reUc or halo) radio emission 
by 1-10 GeV electrons (characterized by a power-law spectral 
distribution Ne{E) cx E~'' with > 3) in the expected OipG) 
magnetic fields. The radio emission appears to correlate with 
merger activity and SXR luminosity of the host cluster. Re- 
cently, excess fluxes beyond the expected thermal emission 
have been discovered in the extreme ultraviolet (EUV; 0.07-3 
keV) and hard X-ray (HXR; 10-100 keV) ranges in several 
clusters by EUVE, BeppoSAX, RXTE, Suzaku, and Integral 
(Rephaeli et al. 2008). The extent and strength of these emis- 
sions has been disputed (e.g. analyses of SwiftlBKI observa- 
tions by AjeUo et al. 2009, and Suzaku observations by Wik et 
al. 2009) and characterization of these radiations as thermal or 
NT remains controversial. However, the undisputed presence 
of radio-emitting leptons indicates that there will be gamma- 
ray emission by these leptons either by inverse Compton (IC) 
scattering of cosmic microwave background (CMB) and other 
soft (IR, optical, SXR) photons, or by NT bremsstrahlimg 
scattering with background charged particles (e.g. Atoyan & 
Volk 2000; Reimer et al. 2004; Blasi et al. 2007). For details, 
see review by Petrosian et al. (2008) and references therein. 

Unlike for the leptons, there is no direct observational ev- 
idence for CR ions in the ICM of any cluster. Previous ob- 
servations by satelUte-based "MeV/GeV" (Sreekumar et al. 
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1996; Reimer et al. 2003) and imaging air-Cherenkov "TeV" 
telescopes (Hattori & Nishijima 2005; Perkins et al. 2006; 
Perkins 2008; Galante et al. 2009; Aharonian et al. 2009a,b; 
Kiuchi et al. 2009; Aleksic et al. 2010), and shower particle 
detectors (Dingus et al. 2005) have provided gamma-ray flux 
upper limits which can constrain populations of hadronic CRs 
from GeV to multi-TeV energies. 

An additional source of GeV gamma rays could be decay 
or annihilation of dark matter particles. Dark matter annihi- 
lation could also contribute to the pool of relativistic leptons. 
However, in view of the currently speculative nature of these 
mechanisms, conventional astrophysical processes should be 
fully understood before attributing observed gamma rays to 
dark matter annihilation. A dedicated analysis of Fenni-LAI 
observations of galaxy clusters in the context of specific dark 
matter models is presented separately (Abdo et al. 2010a). 

In the next section, we describe the selection procedure for 
candidate clusters, ^outlines our analysis of Fenni-LAJi data 
and we present our results, consisting only of upper limits. We 
briefly discuss the constraints these observations place on the 
spectra and distribution of NT particles in the ICM in ^and 
conclude with ^ 

2. CANDIDATE SELECTION 

The selection criteria of clusters with the greatest chance 
of detection at GeV gamma-ray energies are complicated, in- 
volving relativistic ions and electrons, magnetic fields, clus- 
ter merger history, and other properties of the ICM such as 
the thermal gas density and temperature. The gamma-ray 
flux from neutral-pion decay is expected to correlate with 
SXR flux since both emissions are related to the thermal gas 
density (Colafrancesco & Blasi 1998). Therefore, we be- 
gan with the HIFLUCS flux-limited sample of the brightest 
X-ray clusters (Reiprich and Bohringer 2002) and selected 
clusters with the highest ratio of mass-to-distance-squared, 
M /S. To this sample, we added clusters with previously ob- 
served NT emission in other wavebands (e.g. diffuse radio 
features) as an indicator of relatistic lepton populations, and 
also clusters with exceptional X-ray luminosity and tempera- 
ture such as the Bullet cluster (z = 0.296; Clowe et al. 2004), 
MACSJ0717.5H-3745 (z = 0.546; Ebeling et al. 2007), and 
RJ1347.5-1145 (z = 0.451; Schindler et al. 1995). The char- 
acteristics of thirty-three galaxy clusters selected for study in 
LAT data are provided in Table [T] 

3. OBSERVATIONS, ANALYSIS, & RESULTS 

The LAT is a pair-conversion telescope with a precision 
tracker and calorimeter, a segmented anti-coincidence de- 
tector (ACD) which covers the tracker array, and a pro- 
grammable trigger and data acquisition system. The energy 
range of LAT sensitivity spans from 20 MeV to > 300 GeV 
with an angular resolution per single event of approximately 
5.1° at 100 MeV and narrowing to about 0.14° at 10 GeV.^' 
Full details of the instrument, on-board and ground data pro- 
cessing, and other mission-oriented support are given in At- 
wood et al. (2009). 

We searched for high-energy emission from thirty-three 
galaxy clusters using LAT data collected from the commence- 
ment of scientific operations in early-August 2008 to 4 Febru- 
ary 2010. During this period, the LAT operated primarily in a 

Angular resolution is defined here as the 68% containment radius of the 
LAT point spread function averaged over the instrument acceptance and in- 
cluding photons which convert in either the thick or thin layers of the tracker 
array. 



scanning mode (the 'sky survey' mode) that covered the full 
sky every two orbits (i.e., ~ 3 hrs). For operational reasons, 
the standard rocking angle (defined as the angle between the 
zenith and center of the LAT field of view) for survey mode 
was increased to 50° on 3 September 2009 and the data se- 
lection criteria for the present work have been adjusted ac- 
cordingly relative to the IFGL catalog analysis (Abdo et al. 
2010b). We began by selecting all gamma rays of energy 
0.2-100 GeV within a 10° radius around the direction of each 
galaxy cluster in our sample. Only events satisfying the stan- 
dard low -background 'diffuse' class^^ selection criteria corre- 
sponding to the post-launch P6V3 instrument response func- 
tions are accepted into the following analysis. In order to re- 
duce the effects of the so-called 'albedo' gamma rays (from 
interaction of CRs with the upper atmosphere), we remove 
photons arriving from zenith angles > 100° and exclude time 
periods when the rocking angle exceeded 52°. 

The data were prepared using the LAT Science Tools soft- 
ware package and analyzed in the context of diffuse gamma- 
ray emissions and discrete sources previously detected by the 
LAT^^* Diffuse gamma-ray emission from the Milky Way is 
estimated using model templates from gll_iem_v02 . fit. 
The isotropic diffuse component, attributed to both diffuse ex- 
tragalactic gamma-ray emission and residual background of 
charged particles triggering the LAT, has been treated with 
isotropic_iem_v02.txt. In addition, all individual 
objects detected by the Fermi-LAJi appearing in the IFGL cat- 
alog (Abdo et al. 2010b) within a 10° radius of each candidate 
position are included in the sky model as separately-fit point 
sources with fixed position. 

We created detection significance maps of the regions of the 
sky around each cluster using a maximum likelihood analysis 
tool (sourcelike). With the exception of the Perseus and 
Virgo clusters, no significant gamma-ray signal was detected 
by the LAT towards any galaxy cluster in our sample. The 
flux from the Perseus cluster is dominated by variable emis- 
sion from NGC 1275 (probably related to the radio core 3C84; 
Abdo et al. 2009a) while the radio galaxy M87 appears to be 
the primary gamma-ray source of the Virgo cluster (Abdo et 
al. 2009b). 

We determined flux upper limits for each cluster using the 
Neyman construction suitable for experiments with low sig- 
nal counts in the presence of background. Confidence in- 
tervals were created using the Roe-Woodroofe (1999) mod- 
ification to the unified approach proposed by Feldman and 
Cousins (1998). Within the full 10° radius region of interest 
centered on each cluster, we define an inner circular aperture 
for the purpose of comparing the observed counts to the num- 
ber of photons expected from backgrounds alone. The radii of 
the inner apertures are optimized to maximize the signal-to- 
noise ratio given the particular background distributions and 
assumed extent of high-energy cluster emission. The radii of 
photon counting apertures are adjusted to maintain uniform 
signal containment over the 0.2-100 GeV energy range. 

Several clusters in our sample are sufficiently extended to 
be potentially resolved by the LAT, depending on the (as 
yet unknown) spatial distributions of their gamma-ray emit- 



See |http : / /f ermi . gsf c ■ nasa ■ gov/ ssc/ data/ 

lanalys is /documentat ion/Cicerone/ Cicerone_Data/LAT_ 
pp.htmll 

°^ Information regarding the LAT Science Tools package, diffuse mod- 
els, instrument response functions, and public data access is available from 
the Fermi Science Support Center I http : / /f ermi . gsf c . nasa . gov/ 1 
ssc/ I 
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ting regions (see e.g. Pinzke & Pfrommer (2010) and Don- 
nert et al. (2010) for gamma-ray surface intensity distribu- 
tions predicted from numerical simulations). Since soft X-ray 
(bremsstrahlung) emissivity is proportional to the square of 
the thermal gas number density, then assuming that CRs re- 
flect the thermal gas distribution, the expected surface bright- 
ness from hadronic gamma-ray emission should follow the 
X-ray emission. Therefore, we apply a King (/?) surface in- 
tensity profile with core radius rc. 



E(r) = S(0) 



1 + 



2-1 -3/3+0.5 



(1) 



where X-ray profiles are available and can be sufficiently 
well-fit by this form (Chen et al. 2007; Matsushita et al. 2002 
for the Virgo cluster). In addition, we have treated the Coma, 
Fornax, and Virgo clusters with two-dimensional Gaussian 
surface intensity models to assess the stability of the upper 
limits given the uncertain spatial distribution of high-energy 
cluster emission. Tables |2] and [3] provide 95% confidence- 
level flux upper limits in the energy range 0.2-100 GeV as- 
suming unresolved (point-like) or spatially extended hypothe- 
ses, respectively. The gamma-ray flux upper limits presented 
here do not depend strongly on the particular surface intensity 
distributions of cluster emission. 

Integral flux upper limits over a broad energy range depend 
partly upon the assumed spectrum of radiation because the 
LAT effective area increases between 0.2 and 1 GeV. Here, we 
assume a power-law spectrum of high-energy cluster emission 
with photon index =2. The degree to which upper limits 
computed with sourcelike depend upon the underlying 
spectrum of radiation can be estimated by rescaling the obser- 
vation exposures according to the assumed source spectrum. 
For 1.5 < < 3.0, the relative change in photon flux upper 
Umits in the energy ranges 0.2-100 GeV, 0.2-1 GeV, 1-10 
GeV, and 10-100 GeV are 37%, 16%, 2%, and 1%, respec- 
tively, with the LAT being more sensitive to hard-spectrum 
sources. The 1-10 GeV energy band is typically the most 
sensitive in terms of average differential energy flux. 

The flux upper limits to GeV emission from clusters of 
galaxies derived from LAT data are the most stringent to date. 
Upper limits set by EGRET (Reimer et al. 2003) and LAT 
observations are compared to recent model predictions in Fig- 
ure [1] Although the LAT accumulated nearly uniform expo- 
sure over the full sky during the 18-month observation period, 
sensitivity to several candidates including the Ophiuchus and 
Perseus clusters is adversely affected by foreground emission 
from the Galactic plane and other bright gamma-ray sources. 

4. INTERPRETATION 

We now use the flux upper limits provided by LAT ob- 
servations to constrain NT populations of leptons (e*) and 
hadrons (protons and other ions) in the ICM. As described in 
^ gamma rays can be produced from decay of vr^s produced 
in proton-proton interactions between CR ions and ambient 
thermal gas. These collisions inevitably produce secondary 
pairs which contribute to the population of energetic lep- 
tons in the ICM. 

4. 1 . Leptonic Emission Processes 

Both primary and secondary electrons and positrons pro- 
duce radio waves via synchrotron radiation in a magnetic 
field B ~ /LtG and gamma rays either by bremsstrahlung or 



Be- 



by IC scattering of ICM soft photons (energy density Mjoft)- 
For a power-law electron spectrum, A^(7) = C7~'', the ratio of 
gamma-ray flux to radio flux (photon energy = hv) from 

these processes is roughly e^F{e~^) ™ / e^Fie^) <yin/B~ 

and e^F{e^f^ / e^F(e-y)™'' cx Mph/B^ . Thus, for clusters with 
measured radio flux and knowledge of particle and soft pho- 
ton densities, upper limits on gamma-ray fluxes can con- 
strain the volume-averaged value of the magnetic field. We 
use the Coma cluster to illustrate such an analysis. The ra- 
dio spectrum can be fitted by a power law, e^F{e^)"^'^ ~ 

lO"'"* (cMi) ^ ergcm-^s-i for 0.03 GHz <v<\ GHz, requir- 
ing an electron distribution with index /:> ^ 3.5 in the Lorentz 

^ -I 

factor range of 2 x 10^ {j^^ ' < 7 < 5 x 10^ 

yond this energy, the spectrum steepens either as a power 
law with the index higher by one (Rephaeli 1979) or with 
exponential cut-off (Schlickeiser et al. 1987), requiring cor- 
responding steepening of the electron spectrum (Petrosian 
2001). 

For a power-law distribution, the IC flux is given by 

O-p) (£-3) 

e^F(e^) oc e-y- u^ob^^^^^ Mp)^ where A(p) depends only on 
the electron spectral index, thereby favoring IC scattering of 
CMB photons with the highest energy density of mcmb ^ 0.26 
eV cm and ccmb ^ 0.001 eV. However, for p > 3, the con- 
tribution of higher-energy seed photons with lower energy 
density becomes increasingly important. For example, opti- 
cal photons of Esoft = 3 eV and Mopt ^ 0.05 eV cm"-' will have 
equal contribution to CMB for p = 4. Moreover, since the 
electron spectrum steepens for 7 > 10"*, the contribution from 
CMB photons will come primarily from electrons above this 
energy so that (for both the power law or exponential cut-off) 
the main contribution in the GeV range will come from scat- 
tering of optical photons. Similarly, the gamma-ray flux due 
to IR photons will be a factor >100 smaller A detailed 
calculation of the expected IC flux when combined with the 
upper limit set by Fermi implies a lower limit of ^ 0.15/iG 
on the volume-averaged magnetic field of the Coma cluster. 
Furthermore, the non-detection of diffuse gamma-ray emis- 
sion towards any cluster during the first year of LAT observa- 
tions disfavors lepton acceleration efficiencies in intracluster 
shocks >0.001 (Gabici & Blasi 2004). 

The bremsstrahlung flux of the above electrons will pro- 
duce a steeper gamma-ray spectrum e^F(e^) cx e^"'' and will 
require slightly smaller magnetic field lower limit. 

4.2. Hadronic Emission Processes 

For clusters with previously measured gas density (n) and 
temperature (T) profiles, gamma-ray flux upper limits can di- 
rectly constrain the energy density stored in hadronic CR rela- 
tive to the thermal energy density of the ambient ICM gas (e.g. 
Enfilin et al. 1997; Pfrommer & EnBlin 2004), expressed as 



/^CR\_ jNp(E)EdE 



InkT 



(2) 



when averaging over the cluster volume. We assume a uni- 
form power-law spectra of CR hadrons throughout the ICM, 

Klein-Nishina suppression in this case will steepen the gamma-ray spec- 
tram above 20 GeV (electron energy > 10 GeV). IC scattering of more nu- 
merous SXR photons occurs in the Klein-Nishina regime for of energy 
> 100 MeV which will suppress gamma-ray flux above 30 MeV. 
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Np{E) cx £■""'', and a spatial distribution which follows the 
thermal gas. Note that gamma-ray emissivity could be en- 
hanced by a relative over- abundance of CR ions in central re- 
gions where thermal gas density is greatest. The CR-hadron- 
to-thermal energy density ratios for clusters with available X- 
ray data are presented in Table [3] for two possible CR spec- 
tral indices, ap. Using the Coma cluster as an example, we 
find {€cr/€th) < 0.05 compared to the upper limit of < 0.1- 
0.3 obtained from analysis of EGRET observations (Reimer et 
al. 2004). The constraints presented here complement recent 
observations by the H.E.S.S. and MAGIC telescopes which 
have yielded upper limits to {ecR/^Tn) of < 0.2 for the Coma 
cluster (Aharonian et al. 2009b) and < 0.05 for the Perseus 
cluster (Aleksic et al. 2010; simplified model) considering 
populations of multi-TeV CR hadrons in the ICM. CR pop- 
ulations are unlikely to contribute significant pressure which 
would bias X-ray mass estimates of clusters. 

The contraints on hadronic CR populations derived from 
LAT data are in agreement with limits placed by indirect 
methods (Brunetti et al. 2007; Churazov et al. 2008) and 
with the predictions of theoretical models and numerical sim- 
ulations pointing out morphological and spectral difficulties 
(namely, observed radio spectra cut-offs) of explaining large- 
scale radio halos with purely secondary emission (e.g. Blasi 
et al. 2007 and references therein; Donnert el al. 2010). For 
the clusters examined thus far, multiwavelength evidence sug- 
gests that secondary electrons play a minor role in NT emis- 
sion. 

5. CONCLUSIONS 

We have presented gamma-ray flux upper limits for thirty- 
three clusters of galaxies obtained from 18 months of Fermi- 
LAT observations which are robust with respect to uncer- 
tainty in both the spatial extent and spectrum of high-energy 
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TABLE 1 

Summary of the thirty-three galaxy cluster candidates 



Cluster 


; 


b 


z 




Score 




Diffuse radio 


Lx (0. 1-2.4 keV) 


Tx 




(deg) 


(deg) 




(deg) 


(deg) 


(lO'Mg/Mpc^) 




(lO'^ergs-') 


(IceV) 


X-ray flux selection 


3C129 


160.43 


0.14 


0.0223 


0.67 


0.14 


29.1 




2.27 


5.57 


A0754 


239.25 


24.75 


0.0528 


0.40 


0.05 


12.8 




3.97 


9.00 


A 1367 


234.80 


73.03 


0.0216 


0.77 


0.18 


42.7 




1.20 


3.55 


A2199 


62.94 


43.69 


0.0302 


0.46 


0.05 


12.5 




4.20 


4.28 


A2256 


111.10 


31.74 


0.0601 


0.33 


0.10 


8.5 


Halo, Relic (1, 2) 


9.24 


6.83 


A2319 


75.67 


13.58 


0,0564 


0.37 


0.05 


10.9 


Halo (1, 2) 


16,37 


8.84 


A3376 


246.52 


-26.29 


0,0455 


0.36 


0.17 


8.5 




2,16 


4.43 


A3571 


316.32 


28.55 


0,0397 


0.45 


0,05 


14,5 




8,08 


6.80 


Antlia (S636) 


272.94 


19.19 


0,0116 


0.85 


0.29 


31,6 




0,38 


2.06 


AWM7 


146.35 


-15.62 


0,0172 


0.85 


0,10 


45,0 




2,10 


3.70 


Centaurus (A3526) 


302.41 


21,56 


0,0499 


1.24 


0,04 


87,9 




1,19 


3.69 


Coma (A1656) 


58.09 


87,96 


0,0232 


0.80 


0,15 


49,6 


Halo, Relic (1) 


8,09 


8.07 


Fornax (S373) 


236.72 


-53,64 


0,0046 


2.01 


0,36 


168,1 




0,08 


1.56 


Hydra (A1060) 


269.63 


26,51 


0,0114 


1.02 


0,08 


52,5 




0,56 


3.15 


M49 


286.92 


70,17 


0,0044 


1.68 


0,02 


95,5 




0,02 


1.33 


NGC4636 


297.75 


65,47 


0,0037 


1.27 


0,02 


36,3 




0,02 


0.66 


NGC5044 


311.23 


46,10 


0,0090 


0.74 


0,01 


16,6 




0,18 


1.22 


NGC5813 


359.18 


49,85 


0,0064 


1.00 


0,04 


28,9 




0,02 


0.76 


NGC5846 


0.43 


48.80 


0.0061 


0.78 


0.01 


13.3 




0,01 


0.64 


Norma (A3627) 


325.33 


-7.26 


0.0163 


0.89 


0.18 


50.2 




3.59 


5.62 


Ophiuchus 


0.56 


9.27 


0.0280 


0.10 


O.IO 


131.6 


Halo (3) 


12.14 


10.25 


Perseus (A0426) 


150.58 


-13.26 


0.0183 


0.85 


0.03 


49.0 




16.39 


6.42 


Triangulum 


324.48 


-11.63 


0.0510 


0.42 


0.06 


14.7 




12.43 


9.06 


Non-thermal selection 


A0085 


115.05 


-72.06 


0.0556 


0.31 


0.02 




Relic (1,4) 


9.67 


6.51 


A1914 


67.20 


67.46 


0.1712 


0.13 


0.02 




Halo (1,2) 


17.04 


8.41 


A2029 


6.51 


50.55 


0.0767 


0.25 


0.01 




Halo (3) 


17.07 


7.93 


A2142 


44.21 


48.70 


0.0899 


0.24 


0.02 




Halo (4) 


21.05 


8.46 


A2163 


6.75 


30.52 


0.2010 


0.12 


0.03 




Halo(l) 


32.16 


10.55 


A2744 


8.90 


-81.24 


0.3080 








Halo (1) 






Bullet (IE 0657-56) (a) 


266.03 


-21.25 


0.296 








Halo (5) 




14 


MACSJ0717.5+3745 (b) 


61.89 


34,02 


0.546 








Relic (6) 


24.6 


11,6 


Other selection 


RXJ1347.5-1145(c) 


324.04 


48,80 


0.451 










62.0 




Virgo (M87 sub-clump) (d) 


283.78 


74,49 


0,0036 




0,05 











X-ray measurements from Chen et al. 2007 unless stated otherwise. Other X-ray references: (a) Clowe et al. 2004; (b) Ebeling 
et al. 2007; (c) Schindler et al. 1995; (d) Matsushita et al. 2002. Radio references: (1) Giovannini et al. 1999 and references 
therein; (2) Kempner & Sarazin 2001; (3) Govoni et al. 2009; (4) Giovannini & Feretti 2000; (5) Liang et al. 2000; (6) Edge et 
al. 2003. A ACDM cosmology with Hubble constant Hq = 50 km s"' Mpc"' is used for comparison to X-ray data sets. 
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TABLE 2 

95% CONHDENCE-LEVEL GAMMA-RAY FLUX UPPER LIMITS FROM CLUSTERS OF GALAXIES: EMISSION UNRESOLVED BY THE LAT (POINT SOURCE) 



Cluster 


fsrc 


N 














0.2-100 


0.2-100 


0.2-1 


1-10 


10-100 


3C129 


0.54 


3.44 


1.10 


2.04 


1.06 


1.84 


A0754 


0.51 


3.42 


1.10 


1.54 


0.72 


2.15 


A1367 


0.49 


1.65 


0.53 


0.97 


0.93 


2.13 


A1689 


0.50 


4.96 


1.59 


2.14 


0.66 


2.16 


A1914 


0.54 


1.19 


0.38 


0.72 


0.57 


1.86 


A2029 


0.49 


3.28 


1.05 


1.75 


1.20 


2.16 


A2142 


0.52 


2.82 


0.90 


1.59 


0.45 


1.95 


A2163 


0.54 


5.51 


1.77 


2.25 


2.04 


1.99 


A2199 


0.52 


1.12 


0.72 




0.56 


1.91 


A2256 


0.48 


1.96 


0.63 


0.85 


0.70 


1.75 


A2319 


0.51 


0.75 


0.24 


0.42 


0.54 


1.91 


A2744 


0.54 


2.49 


0.80 


1.39 


0.40 


2.01 


A3266 


0.56 


9.23 


2.96 


4.08 


0.61 


1.85 


A3376 


0.57 


9.93 


3.18 


5.18 


1.18 


1.89 


A3571 


0.51 


2.80 


0.90 


1.66 


0.84 


2.13 


A3888 


0.44 


2.41 


0.77 


1.36 


0.72 


2.47 


A85 


0.46 


2.13 


0.68 


1.12 


0.53 


2.35 


AWM7 


0.55 


3.84 


1.23 


1.96 


0.82 


1.82 


Antilia 


0.51 


4.05 


1.30 


2.45 


0.49 


2.13 


Bullet 


0.47 


2.75 


0.88 


1.46 


0.78 


2.24 


Centaurus 


0.51 


8.01 


2.57 


4.12 


1.35 


2.13 


Coma 


0.51 


4.58 


1.47 


1.98 


0.77 


2.01 


Fornax 


0.54 


5.06 


1.62 


3.13 


0.29 


2.00 


Hydra 


0.60 


2.21 


0.71 


0.96 


0.73 


1.81 


M49 


0.52 


2.02 


0.65 


1.13 


0.40 


2.05 


MACSJ0717 


0.54 


6.63 


2.13 


2.89 


1.11 


1.84 


NGC4636 


0.47 


2.96 


0.95 


1.81 


0.45 


2.29 


NGC5044 


0.50 


1.90 


0.61 


0.88 


0.62 


2.18 


NGC5813 


0.52 


10.51 


3.37 


4.22 


1.06 


2.06 


NGC5846 


0.55 


13.02 


4.18 


5.38 


0.69 


1.94 


Norma 


0.47 


1.31 


0.42 


0.82 


0.77 


3.47 


Ophiuchus 


0.54 


26.22 


8.41 


13.32 


2,11 


2.00 


Perseus 


0.68 


89.19 


28.60 


30.34 


28.76 


12.92 


RXJ1347 


0.55 


2.57 


0.82 


1.17 


0.49 


1.98 


Triangulum 


0.54 


2.68 


0.86 


0.90 


0.97 


1.93 


Virgo 


0.62 


14.13 


4.53 


4.97 


3.89 


2.67 



fsrc is the estimated signal fraction captured within the photon counting aperture. and {c~.F{€~^)) correspond to the integral 
photon flux (10~^ ph cm~^ s"') and average differential energy flux (10~^^ erg cm"^ s~'), respectively, within the energy range 
provided (GeV) assuming a power-law spectrum of gamma-ray emission with photon index a^=2. Only photons of energy 
0.4-100 GeV were considered in the analysis of A2199 in order to reduce background from nearby bright sources. 
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TABLE 3 

95% CONFIDENCE-LEVEL GAMMA-RAY FLUX UPPER LIMITS FROM CLUSTERS OF GALAXIES: SPATIALLY EXTENDED EMISSION MODELS 



Cluster 


Spatial model (i?68) 
(deg) 




N 

0.2-100 


0.2-100 


(e^F(£^)) 
0.2-1 


(£^F(e^)) 
1-10 


(€^F(e^)) 
10-100 


dp =2.1 


Qp=2.4 


3C129 


King 


0.56 


3.93 


1.26 


1.84 


1.01 + 


2.61 


0.16 


0.12 


A0085 


King 


0.46 


2.23 


0.71 


1.17 


0.53+ 


2.35 


0.06 


0.04 


A0754 


King 


0.54 


3.31 


1.06 


1.65 


1.06+ 


2.03 


0.35 


0.27 


A1367 


King 


0.48 


1.83 


0.59 


1.03 + 


1.09 


2.18 


0.26 


0.16 


A1914 


King 


0.54 


1.18 


0.38+ 


0.72 


0.57 


1.86 


0.38 


0.24 


A2029 


King 


0.49 


3.18 


1.02 


1.72 


1.19+ 


2.16 


0.11 


0.09 


A2142 


King 


0.52 


2.75 


0.88 


1.59 


0.49+ 


1.95 


0.07 


0.05 


A2163 


King 


0.54 


5.50 


1.76+ 


2.32 


2.17 


1.99+ 


0.81 


0.61 


A2199 


King 


0.51 


1.18 


0.76 




0.54+ 


1.95 


0.14 


0.11 


A2256 


King 


0.47 


1.83 


0.59 


0.81 


0.66+ 


1.79 


0.16 


0.12 


A2319 


King 


0.54 


0.73 


0.23 + 


0.54 


0.58 


2.82 


0.03 


0.02 


A3376 


King 


0.57 


9.69 


3.11 


5.18 


1.08+ 


1.89 


1.21 


0.95 


A3571 


King 


0.50 


3.26 


1.04 


1.85 


0.89+ 


2.17 


0.05 


0.04 


Antlia 


King 


0.52 


4.84 


1.55 


2.86 


0.66+ 


2.09 


1.52 


1.19 


AWM7 


King 


0.55 


3.95 


1.27 


1.92 


0.97+ 


1.82 


0.10 


0.08 


Centaurus 


King 


0.51 


8.15 


2.61 


3.90 


1.33+ 


2.13 


0.09 


0.07 


Coma 


Gauss (0.2) 


0.53 


4.84 


1.55 


2.28 


0.76 


3.03 






Coma 


Gauss (0.4) 


0.52 


4.86 


1.56 


2.36 


0.82 


3.08 






Coma 


Gauss (0.6) 


0.58 


5.12 


1.64 


2.38 


0.82 


4.97 






Coma 


Gauss (0.8) 


0.56 


4.93 


1.58 


2.73 


0.74 


4.88 






Coma 


King 


0.55 


5.14 


1.65 


2.18 


1.11 + 


2.92 


0.05 


0.04 


Fornax 


Gauss (0.2) 


0.51 


4.77 


1.53 


2.61 


0.31 


2.12 






Fornax 


Gauss (0.4) 


0.62 


5.40 


1.73 


2.73 


0.39 


2.72 






Fornax 


Gauss (0.6) 


0.59 


5.73 


1.84 


3.02 


0.43 


2.68 






Fornax 


Gauss (0.8) 


0.62 


5.39 


1.73 


2.61 


0.89 


2.28 






Fornax 


Gauss (1.0) 


0.60 


5.03 


1.61 


2.87 


0.90 


2.30 






Fornax 


King 


0.60 


5.64 


1.81 


2.80 


0.40+ 


2.81 


0.75 


0.59 


Hydra 


King 


0.60 


2.24 


0.72+ 


0.94 


0.85 + 


2.83 


0.28 


0.21 


M49 


King 


0.52 


2.08 


0.67 


1.14 


0.39+ 


2.05 


5.09 


3.98 


NGC4636 


King 


0.46 


2.67 


0.86 


1.28 


0.42+ 


2.34 


3.89 


3.04 


NGC5044 


King 


0.50 


1.87 


0.60 


0.81 


0.62+ 


2.18 


1.58 


1.24 


NGC5813 


King 


0.52 


10.57 


3.39 


4.25 


0.99+ 


2.06 


25.59 


20.03 


NGC5846 


King 


0.55 


13.01 


4.17 


5.38 


0.69+ 


1.94 


13.82 


10.82 


Norma 


King 


0.54 


1.21 


0.39+ 


0.94 


0.83 


3.84 


0.03 


0.02 


Ophiuchus 


King 


0.54 


26.22 


8.41 


14.18 


2.02+ 


1.95 


0.05 


0.04 


Perseus 


King 


0.70 


87.36 


28.01 


28.11 


29.02+ 


17.95 + 


0.27 


0.32 


Triangulum 


King 


0.54 


2.39 


0.77+ 


0.91 


1.03+ 


1.93 


0.07 


0.05 


Virgo 


Gauss (0.2) 


0.62 


14.49 


4.65 


4.93 


4.13 


4.35 






Virgo 


Gauss (0.4) 


0.61 


15.27 


4.90 


5.26 


4.65 


5.13 






Virgo 


Gauss (0.6) 


0.64 


14.97 


4.80 


5.62 


4.72 


4.48 






Virgo 


Gauss (0.8) 


0.64 


15.76 


5.05 


5.62 


4.72 


4.48 






Virgo 


Gauss (1.0) 


0.66 


16.01 


5.13 


5.71 


4.59 


4.35 






Virgo 


Gauss (1.2) 


0.64 


17.03 


5.46 


5.62 


4.72 


4.48 






Virgo 


King 


0.61 


14.89 


4.77 


5.24 


4.33+ 


5.26 


0.17 


0.13 



High-energy cluster emission has been treated with either a two-dimensional Gaussian or King profile spatial model. Gaussian 
models are parameterized by the 68% surface intensity containment radius (Res) and King profiles are fitted to the X-ray surface 
intensity (see text). {ecR/erh) represents the volume-averaged CR-hadron-to-thermal energy density ratio assuming that hadronic 
CRs trace the ambient thermal gas distribution and are described by a power-law spectrum, Np(E) cx Upper limits to 

{^cr/^tii) are computed using the most constraining energy bands (indicated by f). Notes from Table|2]also apply here. 
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Fig. 1. — Photon flux upper limits derived from Feniii'-LAT observations of galaxy clusters (assuming unresolved gamma-ray emission) are compared to 
EGRET results (Reimer et al. 2003) and to recent predictions based on the numerical simulations of Pinzke & Pfrommer (2010) (flux from within one virial 
radius) and Donnert et al. (2010). The LAT integral fluxes presented in Table|2]have been extrapolated to > 0.1 GeV for ease of comparison. 



